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LIGHT  ABSORPTION  BY  ION  ACX)USTIC  TURBULENCE 
IN  LASER-PRODUCED  PLASMAS 


E 


F 


I INTRODUCTION 

The  problem  of  how  laser  light  is  absorbed  in  a laser  produced  plasma  is  of  crucial  impor- 
tance to  the  laser  fusion  program,  as  well  as  being  a very  interesting  scientific  problem  in  its 
own  right.  This  paper  studies  anomalous  absorption  of  laser  light  by  ion  acoustic  turbulence  in 
the  underdense  plasma.  It  is  a continuation  of  our  earlier  work  on  this  subject*,  which  studied 
laser  light  absorption  in  an  unmagnetized  plasma  by  self  generated  ion  turbulence.^  In  this  pa- 
per we  discuss  this  mechanism  further,  and  also  calculate  absorption  by  ion  acoustic  turbulence 
in  a magnetized  plasma.  ^ 

i 

A crucial  characteristic  of  this  absorption  mechanism  is  that  it  occurs  in  the  underdense 
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plasma  where  gradient  scale  lengths  are  expected  to  be  fairly  long,  rather  than  at  the  critical  | 

1 

density  where  the  gradient  scale  length  is  expected  to  be  short.  It  is  interesting  to  note  that 
resonant  absorption,  which  occurs  at  the  critical  surface,  works  best  if  the  density  gradient 
scale  lengths  are  short.  Therefore,  since  these  two  absorption  mechanisms  occur  at  different 
spatial  positions,  there  is  no  reason  why  they  cannot  both  occur  simultaneously.  A schematic 

of  the  expected  plasma  density  profile,  with  the  regions  of  resonant  absorption,  and  absorption  j 
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by  ion  acoustic  turbulence  is  shown  in  Fig.  1.  The  units  are  free  space  wavelength  of  the  laser 
light.  * 

Another  interesting  aspect  of  absorption  by  ion  acoustic  turbulence  is  that  a single  fluid 
theory  is  valid  for  all  electrons.  Although  the  classical  mean  free  path  of  an  energetic  electron 
is  long  compared  to  a gradient  scale  length,  this  mean  free  path  is  anomalously  reduced  in  the 
Manuicript  tubmitted  Auguit  6,  1977. 
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presence  of  turbulence  and/or  a magnetic  field. Thus  it  is  possible  that  a single  fluid  theory, 
with  anomalous  transport,  but  without  special  treatment  for  energetic  electrons,  can  give  an 
adequate  description  of  electron  thermal  transport  in  a laser  produced  plasma. 

Section  II  of  this  paper  reviews  recent  data  on  laser  light  absorption  at  a number  of 
different  laboratories  and  shows  that  resonant  absorption  cannot  account  for  more  than  about 
one  fourth  to  one  half  of  the  total  absorption.  Section  III  discusses  anomalous  absorption  in  an 
unmagnetized  plasma,  while  Sections  IV  and  V discuss  absorption  in  a magnetized  plasma. 
Section  VI  shows  how  our  results  scale  with  laser  wavelength,  and  Section  VII  discusses  how 
to  solve  for  the  dynamics  of  the  problem  instead  of  the  steady  state. 

II  THE  NEED  FOR  ANOTHER  ABSORPTION  MECHANISM 

Before  detailing  our  proposed  absorption  mechanism,  we  will  briefly  review  some  laser 
plasma  experiments  and  then  show  that  the  three  absorption  mechanisms  currently  considered, 
inverse  Bremsstrahlung,  resonant  absorption  and  parametric  instabilities  are  not  sufficient  to 
explain  the  measured  laser  light  absorption.  There  are  now  many  different  experiments,  from 
many  different  laboratories  which  show  fractional  light  absorption  of  order  50%  or  more.^ 
Several  of  these  measurements  show  that  the  fractional  absorption  is  not  strongly  dependent 
on  either  tilt  angle  of  the  target'*”  or  polarization  of  the  laser  light.'*  Some  of  these, and 
also  other  experiments'^’’^  measure  the  density  gradient  scale  length  in  the  underdense  re- 
gion. For  neodymium  laser  experiments,  underdense  scale  lengths  are  typically  between  10  and 
60  microns, ^ where  here  the  underdense  scale  length  is  defined  as  the  distance  between 
the  critical  and  quarter  critical  surfaces.  For  COj  laser  experiments,  they  are  typically  larger 
than  50  microns.  ” 

Experiments  at  the  Naval  Research  Laboratory  have  consistently  shown  high  light  ab- 
sorption on  experiments  without  prepulses.  For  instance  Ref  4 shows  absorption  between  fifty 
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and  seventy  percent  which  is  insensitive  to  polarization  and  target  tilt  angle  for  tilt  angles  less 
than  about  50°-60®.  Furthermore,  experiments  in  the  Soviet  Union  with  a neodymium  laser 
have  also  shown  fifty  to  seventy  percent  absorption. 

Experiments  with  CO2  laser  produced  plasmas  show  high  absorption  too.  For  instance 
the  experiments  of  Donaldson  et  al^  show  absorption  of  around  90%.  The  experiments  of 
Richardson  show  absorption  of  about  SO  percent  when  the  target  is  in  the  focal  plane  of  the 
lens.  Furthermore,  this  absorption  is  relatively  independent  of  the  target  tilt  angle  for  tilt  an- 
gles of  less  than  about  thirty  degrees.  Finally,  the  experiments  of  Giovanielli  *et  al^*’ also  show 
that  the  fractional  absorption  is  roughly  45%. 

One  very  interesting  series  of  experiments  is  that  reported  by  Van  Kessel  et  al.^  They 
showed  that  for  a neodynium  laser  produced  plasma,  the  absorption  depends  only  weakly  on 
laser  light  irradiance,  but  depends  quite  strongly  on  the  focal  conditions.  If  the  target  is  in  the 
focal  plane  of  the  lens,  the  absorption  maximizes  at  a value  of  between  forty  five  and  sixty  per- 
cent. As  the  target  was  moved  in  front  of  or  in  back  of  the  focal  plane,  the  fractional  light  ab- 
sorption fell  off.  A similar  result  was  found  in  a recent  eperiment  at  the  Naval  Research 
Laboratory.**  There,  the  laser  irradiance  was  varied  between  4x10*^  W/cm^  and 
2xl0'*W/cm^  at  focus.  The  absorption  was  relatively  constant  at  between  fifty  and  sixty 
percent.  However  it  fell  off  when  the  target  was  moved  out  of  the  focal  plane.  For  instance 
two  shots  at  the  same  intensity,  but  different  focal  conditions  had  markedly  different  fractional 
absorption.  Finally,  we  note  that  this  same  dependence  of  fractional  absorption  on  focal  condi- 
tions was  found  for  a CO2  laser  produced  plasma  also. ' * 

These  results  seem  to  be  consistent  with  some  experiments  made  at  the  Lawrence  Liver- 
more Laboratory^  on  absorption  as  a function  of  laser  irradiance.  In  Ref.  8,  a parylene  disk  was 
irradiated  with  a Nd  laser  pulse  of  between  S and  IS  Joules  with  a duration  of  between  SO  to 
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ISO  picoseconds.  The  fractional  absorption  as  a function  of  irradiance  is  shown  in  Fig.  10a  of 
Ref.  8.  At  the  highest  irradiance,  fractional  absorption  is  greater  than  fifty  percent.  Since  the 
irradiance  in  Fig.  10a  of  Ref.  8 varies  by  more  than  two  orders  of  magnitude,  while  the  irradi- 
ance at  focus  can  only  vary  by  at  most  a factor  of  nine,  much  of  the  variation  in  irradiance 
must  correspond  to  a variation  in  target  position  with  respect  to  the  focal  plane.  Thus  a rea- 
sonable conclusion  appears  to  be  that  in  the  Livermore  experiment,^  absorption  at  best  focus 
is  around  fifty  percent,  consistent  with  the  other  experiments  we  have  cited. 

We  will  now  discuss  absorption  mechanisms  currently  invoked  to  explain  these  data. 
The  favorite  mechanism  now  is  probably  resonant  absorption.*^  If  the  laser  light  wave 
vector  is  not  perpendicular  to  the  target,  and  the  electric  field  is  polarized  in  the  plane  of  in- 
cidence (p  polarization),  then  the  electric  field  perpendicular  to  the  plane  of  the  target  can  tun- 
nel from  the  point  of  reflection  through  to  the  critical  density  and  excite  an  electrostatic  plas- 
ma wave.  Since  this  electrostatic  wave  cannot  propagate  out,  this  energy  is  ultimately  deposited 
into  the  plasma. 

Recent  calculations  show  that  the  absorption  as  a function  of  target  tilt  angle  maximizes 
at  about  fifty  percent. If  the  electric  field  is  polarized  perpendicular  to  the  plane  of  in- 
cidence however  (s  polarization),  there  is  no  resonant  absorption.  It  might  be  possible  to  ex- 
plain (with  resonant  absorption)  experiments  which  show  that  fractional  absorption  is  not 
strongly  dependent  on  tilt  angle  or  polarization.  One  posibility  is  that  laser  beams  focused  by 
real  focusing  lenses  subtend  a substantial  range  of  angles  of  incidence  near  focus.  One  might 
also  argue  that  the  critical  surface  could  be  somewhat  wrinkled.  Experiments,^^  two  dimen- 
sional simulations^^  (for  instance  see  the  figures  of  Ref.  25)  and  theory  all  indicate  this 
possibility.  Thus  the  angle  and  polarization  dependence  could  be  partially  averaged  out.  The 
point  to  keep  in  mind  however,  is  that  each  effect  which  averages  out  the  angle  and  polariza- 
tion dependence  will  probably  also  reduce  the  overall  absorption  efficiency.  For  instance,  no 


matter  how  the  critical  surface  is  wrinkled,  at  each  point  on  it,  the  laser  light  has  two  possible 
transverse  polarizations;  one  with  the  electric  vector  parallel  to  the  plane  of  the  surface,  and 
one  with  a component  of  the  electric  vector  perpendicular  to  the  plane  of  the  surface.  Only 
the  latter  polarization  can  give  rise  to  resonant  absorption.  Thus  while  a wrinkled  critical  sur- 
face may  completely  blur  the  distinction  between  sand  p polarization,  the  fractional  absorption 
will  be  cut  in  half.  Furthermore,  averaging  out  the  angular  dependence  will  reduce  the  overall 
absorption  efficiency  further  still.  Thus,  if  some  mechanism  (for  instance  surface  wrinkling), 
averages  out  the  angle  and  polarization  dependence  of  resonant  absorption,  the  absorption 
efficiency  will  probably  be  reduced  from  about  50%  to  about  15%-20%.  Clearly  this  is  too  low, 
by  at  least  a factor  of  two  to  four,  to  explain  the  numerous  experiments  we  have  cited. 

Recently,  experiments  at  Livermore  have  conclusively  shown  that  resonant  absorption 
does  in  fact  play  an  important  role  in  the  total  absorption  picture.^^  With  the  use  of  a box 
calorimeter,  very  accurate  measurements  of  fractional  absorption  as  a function  of  angle  and  po- 
larization could  be  made.  As  shown  in  Fig.  4 of  Ref.  27,  for  s polarization,  the  fractional  ab- 
sorption at  0 (tilt  angle)  — 0°  was  about  30%.  Fractional  absorption  was  quite  independent  of  0, 
it  fell  gradually  to  about  15%  at  # — 60°.  As  shown  in  Fig.  5 of  Ref  27,  for  p polarization,  the 
same  qualitative  features  were  observed  for  roughly  0 < 9 < 20 “and  45°  < 0 < 60°.  Howev- 
er between  about  20°  and  45°  the  fractional  absorption  for  p polarization  has  a maximum.  For 
instance,  at  » — 30°,  the  fractional  absorption  for  p polarization  varies  between  about  35%  and 
45%,  while  for  s polarization  it  varies  from  about  25%  to  slightly  over  30%.  Thus  resonant  ab- 
sorption is  clearly  playing  an  important  role. 

However  even  this  experiment  appears  to  show  that  other  absorption  mechanisms  are 
present.  One  obvious  question  is  what  is  causing  such  30%  absorption  for  s polarized  light.  The 
fact  that  there  is  such  a distinct  maximum  at  0 — 30°  for  p polarized  light  seems  to  argue 
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against  large  scale  averaging  out  of  angles  or  polarizations.  One  could  in  fact  argue  from  this 
data,  that  there  is  some  underlying  absorption  mechanism  which  is  nearly  independent  of  po- 
larization and  angle,  that  is  roughly  the  experimental  curve  for  s polarization.  However  if  con- 
ditions are  right,  the  absorption  can  further  be  enhanced  by  resonant  absorption. 

Thus  our  conclusion,  based  upon  our  examination  of  the  generally  available  data,  is  that 
while  resonant  absorption  does  play  a significant  role  in  the  overall  picture,  there  is  a clear  and 
pressing  theoretical  need  for  some  other  mechanism  also. 

The  other  theoretical  mechanisms  discussed  up  to  now  have  been  inverse  bremsstrahlung 
and  parametric  instabilities.^*  Both  of  these  mechanisms  are  faced  with  severe  and  possi- 
bly catastrophic  difficulties  in  explaining  light  absorption  in  a laser  produced  plasma.  These 
have  been  discussed  at  some  length  in  the  literature  and  we  will  touch  only  briefly  upon 
them. 

Inverse  bremsstrahlung  is  simply  collisional  damping  of  the  laser  light  as  it  propagates 
through  the  underdense  plasma.  For  intense  laser  light,  temperatures  in  the  underdense  plas- 
ma are  so  high  that  the  collisional  damping  is  too  small  to  give  significant  absorption  in  the 
scale  lengths  relevant  to  a laser  fusion  plasma. 

Parametric  instabilities,  either  the  decay  or  oscillating  two  stream  can  only  occur  near 
the  critical  density.  However  it  is  now  well  established  from  theory, numerical 
simulations, and  microwave  plasma  experiments^^  that  deposition  of  laser  momentum  at 
the  critical  surface  gives  rise  to  a shock  like  discontinuity  at  the  critical  surface  whose  width  is 
about  half  of  the  free  space  wavelength  of  the  laser  light.  This  short  scale  length  around  the 
critical  surface  practically  rules  out  parametric  instabilities  as  an  important  absorption 

mechanism  in  a laser  produced  plasma. 
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A simple  derivation  of  the  homogeneous  threshold  for  the  decay  instability  can  be  given. 
We  show  in  fact  that  for  those  short  scale  lengths  near  the  critical  density,  which  are  appropri- 
ate for  a laser  produced  plasma,  the  decay  instability  is  not  expected  to  be  excited.  The  deriva- 
tion follows  closely  that  given  by  Klein  et  al.^^  More  rigorous  derivations  give  the  same  gen- 
eral results.^*  Say  that  the  density  gradient  (whose  scale  ength  is  L)  is  in  the  x direction  and 
the  laser  light  electric  field  is  in  the  y direction.  Then,  electron  plasma  waves  and  ion  acoustic 
waves  with  wave  number  q in  the  y direction  are  excited.  However  they  are  only  excited  in 
the  spatial  region  where  the  ipn  acoustic  wave,  plasma  wave  and  laser  light  are  all  in  phase. 

If  K is  the  wave  vector  of  the  laser  light,  and  subscripts  / and  e denote  ion  acoustic  and  elec- 
tron plasma  wave  respectively,  the  length  of  the  region  of  phase  coherence,^,  is  defined 
roughly  as 

Since  the  density  jump  scale  length  is  shorter  than  a laser  wavelength,  we  neglect  the 
variation  of  :K^  in  Eu.  (1).  The  dispersion  relation  for  the  ion  wave  is  given  by 


where  Cj  is  the  sound  speed.  Thus  there  is  no  variation  in  with  space  either.  The  plasma 
wave  dispersion  is 


o)^  = (x)  + 3 q^g  + qjf  Vl 


(3) 


where  is  the  electron  thermal  velocity.  At  x *=0,  the  plasma  wave  is  assumed  to  have 
q^^  - 0.  Solving  for  as  a function  of  x and  integrating  Eq.  (1),  we  find 


if 


1 V 

3 /1/3  ^ 

2 2/3  {U^  (x  ”0  ) 


2/3 


(4) 


In  order  for  the  inhomogeneous  plasma  to  be  unstable,  the  waves  have  to  grow 


significantly  while  they  are  in  the  interaction  region,  or 

> 1 


(5) 
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where  y is  the  growth  rate,^’ 


1 

00 

where  Vg^  is  the  electron  oscillatory  velocity  amplitude  in  the  field  of  the  laser  ligi”.  In  Eq. 
(5),  T,  (t,)  is  the  time  the  ion  (electron)  wave  spends  in  the  unstable  region.  Since  the  ion 
wave  convects  out  with  the  flow-  velocity  V, 

r,  ^^IV. 

(7) 

The  electron  wave  refracts  out  of  the  unstable  region.  A simple  calculation  shows  that 

(8) 

Combining  Eqs.  (4-8),  we  find  that  the  condition  for  instability  is 

j_^,  > 1 

16  (9a) 

Assuming  k takes  on  its  maximum  value,  Ej  = V and  L ' as  shown  in 

Ref.  33,  the  condition  for  decay  instability  on  the  density  jump  becomes 

y ^ 

0.1  > 1. 

K (9b) 

For  temperature  above  a few  hundred  electron  volts,  this  condition  cannot  be  satisfied  for 
^os  decay  instability  at  least  is  stabilized  by  the  short  gradient  scale  length. 

Thus  of  the  three  absorption  mechanisms  generally  invoked,  resonant  absorption  is  clear- 
ly the  most  viable.  However,  equally  clearly,  experiments  and  theory  show  that  it  cannot  ac- 
count for  more  than  about  one  third  of  the  measured  absorption.  The  problem  then  is  to  find 
the  cause  for  the  remaining  absorption.  One  possible  approach  is  to  see  if  there  is  any  way 
resonant  absorption  can  be  enhanced.  For  instance,  it  has  recently  been  shown  that  in  a mag- 
netized plasma,  resonant  absorption  can  occur  at  W “ 0 by  coupling  to  an  upper  hybrid  wave 
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rather  than  a plasma  wave.^*  In  this  paper,  we  take  another  approach  and  attempt  to  exam- 
ine a completely  different  absorption  mechanism.  It  is  enhanced  inverse  bremsstrahlung  result- 
ing from  self  induced  ion  density  fluctuations.'  Our  conclusion  is  that  this  is  indeed  a viable 
scheme  which  shows  great  promise  for  explaining  the  high  fractional  absorption  found  in 
many  laser  plasma  experiments. 

Ill  ABSORPTION  BY  SELF  INDUCED  ION  FLUCTUATIONS 
IN  AN  UNMAGNETIZED  PLASMA 

In  this  section  we  review  and  discuss  our  earlier  work  on  laser  light  absorption  by  ion 
acoustic  fluctuations  in  an  unmagnetized  plasma.'  The  idea  is  that  as  energy  is  absorbed,  the 
laser  light  energy  flux  is  converted  into  electron  thermal  energy  flux.  In  order  to  maintain 
charge  neutrality,  a return  current  of  low  velocity  electrons  must  flow  toward  the  laser.  This 
return  current  can  drive  unstable  ion  acoustic  waves  which  also  propagate  toward  the  laser. 
These  ion  acoustic  waves  have  two  principal  effects. 

First  of  all,  they  impede  electron  thermal  transport,  by  anomalously  reducing  the  electron 
mean  free  path;  and  secondly,  they  give  rise  to  enhanced  absorption  of  the  laser  light  as  it  en- 
counters the  ion  density  fluctuations.'  The  theory  of  both  of  these  processes  have  been  given 
in  Refs.  1 and  2 and  we  will  not  detail  it  here. 

The  remainder  of  this  section  will  discuss  the  results  of  these  papers.  In  Ref.  1,  the 
steady  state  fluid  equations  were  integrated  from  x ^0  backwards  toward  the  laser.  Coupled  to 
the  equation  for  the  fluid  quantities  (density,  velocity,  temperature,  and  heat  flux),  were  steady 
state  equations  for  the  laser  light  and  also  for  the  ion  acoustic  fluctuations. 

Many  of  the  fluid  quantities  at  x — 0 could  be  eliminated  by  assuming  that  for  small  posi- 
tive X,  there  is  a density  jump  of  the  type  found  by  Lee  et  al.^^  (Recall  that  the  laser  is  at 
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X — — «>  and  the  integration  is  toward  negative  x).  In  terms  of  the  density  at  x — 0,  the  calcu- 
lations of  Lee  et  al^^  specify  the  fluid  velocity  relative  to  the  ion  acoustic  velocity  and  also 
specify  the  electron  oscillating  velocity  in  the  laser  field  relative  to  the  electron  thermal  veloci- 
ty. Thus  choosing  a lower  shelf  density  and  laser  irradiance  is  equivalent  to  specifying  the 
velocity  and  temperature  at  x “ 0.  The  remaining  quantities  to  determine  were  ion  tempera- 
ture and  electron  energy  flux  Q.  The  electron  energy  flux  was  chosen  so  that  Q(  — oo)  — o, 
and  absorption  was  studied  as  a function  of  ion  temperature  and  density  at  x 0. 

To  solve  for  the  laser  light  absorption,  it  is  necessary  to  know  the  spectrum  of  the  ion 
acoustic  turbulence  not  only  in  wave  number,  but  also  in  angle.  The  reason  is  that  the  return 
current  propagates  anti-parallel  to  the  Poynting  vector  of  the  laser  light.  Thus  the  unstable  ion 
acoustic  waves  propagate  mostly  in  this  direction.  However  only  a density  fluctuation  with  a 
component  of  its  wave  vector  parallel  to  the  laser  light  electric  field  give  rise  to  absorption.^ 
Thus  one  must  either  do  a full  two  dimensional  calculation  or  else  one  must  assume  some  an- 
gular spectrum.  We  chose  the  latter  course  because  earlier  experience^’  has  shown  that  a full 
two  dimensional  calculation  is  extraordinarily  complex,  but  it  adds  very  little  to  the  basic  phy- 
sics. 


The  results  of  our  calculations  were  in  many  ways  very  encouraging.  We  found  high  ab- 
sorption which  was  not  strongly  dependent  on  either  polarization  or  target  tilt  angle.  Also  the 
flux  limits  calculated, / = Q/nmy^  were  around  0.1  which  were  reasonably  consistent  with  ex- 
periments on  layered  targets.*^  Furthermore,  this  absorption  mechanism  does  not  tend  to  pro- 
duce a nonthermal  tail  on  the  electron  distribution  function,  but  rather  tends  to  deposit  energy 
into  particles  whose  velocity  is  |fl/k|  where  fl  is  the  laser  light  frequency  and  k is  the  wave 
number  of  the  density  fluctuation.  This  velocity  turns  out  to  be  about  two  to  four  times  the 
electron  thermal  velocity. 
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We  now  turn  to  some  of  the  difficulties  associated  with  this  absorption  mechanism.  Probably  the 
main  problem  with  the  mechanism,  as  formulated  in  Ref.  l,is  that  it  requires  rather  long  plas- 


I 


ma  lengths  in  the  underdense  region.  For  instance  for  the  calculation  shown  in  Fig.  la  of  Ref. 
1, about  60jli  of  Nd  laser  produced  plasma  was  required  at  half  the  critical  density.  In  a short  pulse 
experiment,  it  is  not  at  all  clear  that  the  plasma  can  expand  so  far.  Even  in  a long  pulse  ex- 


periment, 60/x  of  plasma  at  half  the  critical  density  will  almost  certainly  give  rise  to  strong 

i 

Brillouin  backscatter.^^'^'*  Thus,  before  this  mechanism  can  be  really  an  explanation  for  i 

significant  absorption,  the  necessary  scale  lengths  of  the  underdense  shelf  must  be  reduced. 


There  are  many  ways  this  can  be  accomplished.  In  Ref.  1,  the  ion  acoustic  turbulence 
was  started  out  at  a very  low  level.  As  is  apparent  from  Fig.  la  of  Ref.  1,  about  half  of  the 
necessary  length  was  required  simply  to  build  up  the  turbulence.  If  the  turbulence  started  out 
at  for  instance  ~ 10  at  x =0,  the  length  required  for  absorption  would  be  cut  in 

half. 


There  are  many  reasons  why  e could  be  this  large  at  x =0.  First  of  all,  if  there  is 
any  transport  inhibition  on  the  high  density  side  of  the  critical  surface,  the  ion  acoustic  fluc- 
tuations responsible  for  it  could  propagate  to  the  low  density  side,  initializing  the  turbulent 
fluctuations  there  at  some  large  value.  Secondly,  the  critical  surface  itself  might  be  unstable. 
There  have  been  several  theoretical  and  numerical  investigations  of  instabilities  which  wrinkle 
the  critical  surface. Also,  experiments  at  NRL^^  show  that  the  critical  surface  is  wrin- 
kled with  a mean  angle  of  the  wrinkles  of  about  15°.  This  wrinkling  of  the  critical  surface 
could  certainly  initialize  ion  density  Auctions  of  e4>ITf  ( “Sw/n)  -—0.1.  Thirdly,  many  au- 
thors have  speculated  that  Langmuir  solitons  could  be  excited  at  the  critical  surface  of  a laser 
produced  plasma.'**  '**’  This  also  could  initiate  ion  density  fluctuations  there.  Fin- lly,  it  is 
worth  pointing  out  that  the  solutions  of  Lee  et  al**  show  that  the  low  density  shelf  is  not  a un- 
iform plasma;  rather  it  is  a plasma  with  a large  amplitude,  non  sinusoidal  density  fluctuation 
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whose  wave  length  is  about  one  half  of  the  free  space  laser  light  wavelength.  Surely,  this  large 
amplitude  density  fluctuation  could  also  initiate  ion  acoustic  waves  at  - 0.  Thus  it  is  reason- 
able to  assume  that  the  long  region  in  which  the  turbulence  is  building  up  will  be  eliminated. 
This  reduces  the  necessary  length  by  about  a factor  of  two,  say  from  60jii  to  SOfi 

The  length  required  can  be  still  further  reduced  by  reducing  the  absorption.  For  instance 
if  the  boundary  condition  on  Q were  imposed  at  jc  — — i5#x  rather  than  at  jf  — — <»,  the  prob- 
lem could  be  solved  the  same  way,  but  smaller  values  of  0(x  —0)  (and  therefore  of  fractional 
absorption)  would  be  found.  For  instance  reducing  the  length  by  a factor  of  two,  from  30m  to 
15m  would  probably  just  reduce  the  fractional  absorption  by  a factor  of  roughly  two.  Since  Ref. 
1 found  very  high  absorption,  typically  forty  to  sixty  percent,  the  fractional  absorption  can  be 
reduced  by  a factor  of  two  and  this  will  still  be  a very  important  absorption  mechanism.  Thus 
we  conclude  that  length  scales  can  be  greatly  reduced  from  those  shown  in  Ref.  1 while  ab- 
sorption remains  high. 

Another  possible  problem  is  that  the  mechanism  in  Ref.  1 only  works,  for  the  boundary 
conditions  imposed,  if  TjTj  is  larger  than  about  15.  The  reason  is  that  as  the  ion  temperature 
increases,  enough  ions  have  velocity  = a where  a ^ 1,  comparable  to  the  wave 

phase  velocity  so  that  ion  Landau  damping  is  important  and  the  unstable  region  is 

too  small.  While  temperatures  ratios  as  high  as  15  to  20  are  not  unreasonable,  and  have  in  fact 
been  observed  in  fluid  simulations  of  laser  produced  plasmas,*^  the  absorption  mechanism 
would  be  more  universal  if  this  temperature  ratio  could  be  reduced.  We  find  that  the  tempera- 
ture ratios  for  efficient  absorption  can  in  fact  be  significantly  reduced  if  the  ion  charge  Z > 1. 
The  reason  is  that  the  wave  phase  velocity  increases  to  while  the  ion  velocity 

remains  a so  that  the  importance  of  ion  Landau  damping  is  reduced.  In  a later  sec- 

tion, we  show  that  if  Z - 2.7,  corresponding  roughly  to  a CHj  plasma,  then  there  can  be  ab- 
sorption with  temperature  ratios  as  low  as  five. 
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A final  problem  with  this  absorption  mechanism  is  that  it  depends  on  a number  of 
parameters  which  cannot  be  a-priori  specified.  These  include  the  lower  shelf  density,  the  ion 
temperature,  and  the  boundary  condition  on  0 at  x — — <».  In  a magnetized  plasma,  the  situa- 
tion is  even  worse,  because  the  magnetic  field  must  also  be  specified.  In  a later  section  we  will 
show  how  these  parameters  and  boundary  conditions  can  be  specified  by  solving  for  the 
dynamics  of  the  problem  rather  than  the  steady  state. 

IV  ABSORPTION  IN  A MAGNETIZED  PLASMA, 

QUAUTATTVE  DISCUSSION 

It  is  now  well  established  by  both  theory^®  and  experiment*^  that  magnetic  fields  of 
order  a megagauss  can  be  generated  in  a laser  produced  plasma.  Thus  there  is  excellent 
motivation  for  studying  anomalous  absorption  by  ion  acoustic  waves  in  a magnetized  plasma. 

In  a magnetized  plasma  the  absorption  process  discussed  in  the  previous  section  can  be 
enhanced  by  a very  considerable  amount  as  will  be  discussed  shortly.  Let  us  for  now  assume 
slab  geometry  so  that  the  magnetic  field  is  in  the  z direction,  and  gradients  are  in  the  x direc- 
tion. Then,  it  was  shown  in  Ref.  3 that  ion  acoustic  waves  which  propagate  principally  in  the  y 
direction  are  driven  unstable.  The  growth  rate  is  proportional  to  dTJdx  and  the  electron  ther- 
mal energy  flux  in  the  x direction  is  proportional  to  dTJdx. 

Fluid  equations,  equations  for  the  ion  acoustic  turbulence,  and  equations  for  the  laser 
light  are  integrated  back  from  x - 0 to  x - -<».  For  small  positive  x,  a density  jump  of  the 
type  found  by  Lee  et  al**  is  once  again  assumed.  Thus  the  fluid  velocity  and  temperature  can 
be  determined  in  terms  of  lower  shelf  density  and  laser  irradiance.  We  have  found  that  a 
different  boundary  condition  on  Q(x  - -«>)  must  be  imposed.  The  reason  is  that  if  0 starts 
out  positive  at  x — 0,  it  always  ends  up  negative  at  x — — «»  no  matter  how  parameters  are 
varied.  That  is,  fi  at  x = -<»  can  never  be  zero. 
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Actually  an  outward  energy  flux  for  large  negative  x may  even  be  a more  reasonable  boundary 
condition  than  Q(x  - -oo)--0  as  assumed  in  Ref.  1.  For  instance  if  the  plasma  for  negative  x 
expands  as  a rarefaction  wave, 


n “ rig  exp 


(10a) 


where  Cj  - (Z  T^/M)  , then  the  fluid  energy  flux  to  the  left  at  x - 0 is 

If) MlttgMC]  -SllrigZCJ,  ~ -3ngMC^. 


(10b) 


(11) 


However  one  can  integrate  Eqs.  (10)  to  get  the  rate  of  change  of  energy  in  the  region 
X < 0.  The  energy  in  this  region  increases  faster  than  the  energy  flux  (Eq.  (11))  times  the 
area.  Thus  additional  outward  energy  flux  is  needed  to  feed  the  isothermal  rarefaction  wave. 
The  total  incident  energy  flux  needed  is  -ArigMC^.  Thus  if  the  outward  electron  thermal  en- 
ergy flux  at  X “ -oo  is  —rigMC^,  the  total  outward  energy  flux  will  be  just  sufficient  to  propa- 
gate an  isothermal  rarefaction  wave.  Thus  the  boundary  condition  we  assume  is 


0(x*=-oo)  = — iigMCl 


(12) 


Another  difference  between  the  calculations  presented  in  the  next  section  and  those  in 
Ref.  1 is  that  we  generally  assume  20%  resonant  absorption  at  the  density  jump.  (However  a 
few  calculations  were  done  without  resonant  absorption  as  a comparison).  One  advantage  of 
this  absorption  mechanism  is  that  it  does  not  compete  with  resonant  absorption  but  co-exists 
very  easily  with  it.  This  contrasts  with  the  case  for  parametric  instabilities.  Parametric  insta- 
bilities and  resonant  absorption  can  occur  only  at  the  critical  density.  However  steepening  at 
the  critical  density  tends  to  enhance  resgnant  absorption  and  reduce  the  effect  of  parametric 
instabilities.  Broadening  the  profile  at  the  critical  density  does  the  opposite.  Thus  resonant  ab- 
sorption and  parametric  instabilities  compete  with  each  other. 
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Absorption  by  self  generation  of  ion  turbulence  and  resonant  absorption  do  not  compete 
with  each  other  however,  since  they  occur  at  different  locations  in  the  plasma.  If  the  density 
profile  is  steep  near  the  critical  surface,  but  broad  in  the  underdense  region  as  shown  in  Fig.  1, 
then  both  processes  can  occur  simultaneously. 

There  are  several  reasons  why  this  absorption  mechanism  can  be  enhanced  in  a magnet- 

dTe  . 

ized  plasma.  First  of  all  for  a given  Q(x  -0);  the  smaller  e<l>IT,  is,  the  larger  — is,  and 

therefore  the  larger  the  growth  rate  is.  Thus,  even  if  e<l)ITg  is  assumed  very  small  at  x -0,  the 
growth  rate  is  very  large,  and  the  transient  region  of  small  e<i>ITg  is  effectively  eliminated. 

Secondly,  if  the  laser  light  is  polarized  with  its  electric  field  vector  E in  the  y direction,  it 
and  the  wave  numbers,  k,  of  the  ion  acoustic  turbulence  are  aligned  for  maximum  absorption. 
Thus  there  is  no  need  for  making  any  assumptions  at  all  about  the  angular  spectrum  of  the 
turbulence.  In  fact  a simple  calculation  shows  that  the  same  e^/Tg  and  |k|  gives  rise  to  about  6 
times  as  much  absorption  if  k and  E are  aligned  than  it  does  for  the  angular  spectrum  assumed 
in  Ref.  1.  Thus  here  we  find  roughly  equivalent  absorption,  but  with  much  smaller  values  of 
turbulent  field  strength  and  scale  length.  (Of  course  if  E is  in  the  z direction,  then  one  must 
again  make  assumptions  about  the  angular  width  of  the  spectrum.) 

Thirdly,  the  ion  acoustic  wave  propagates  in  the  y direction  rather  than  the  negative  x 
direction.  Thus  its  group  velocity  in  the  x direction  is  just  the  fluid  speed,  rather  than  the  fluid 
speed  plus  the  sound  speed.  Therefore  the  spatial  growth  rate  yl  is  larger  in  the  magnet- 
ized plasma  case  even  if  y is  unchanged.  To  summarize,  the  effect  of  a magnetic  field  can 
significantly  improve  the  absorption  of  laser  light  by  self  generation  of  ion  acoustic  turbulence. 

We  conclude  this  section  with  some  speculations  on  the  geometry  of  the  magnetic  field  at 
the  critical  surface.  Let  us  take  the  case  of  a focused  laser  beam  striking  a slab  target.  Then 
the  geometry  is  basically  cylindrical;  the  magnetic  field  lines  are  circles  as  shown  in  Fig.  2a. 
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The  polarization  of  the  laser  light  can  be  either  up  and  down  (s  polarization)  or  right  to  left  (p 
ploarization).  It  is  not  clear  whether  the  temperature  gradient  is  principally  axial  or  principally 
radial.  If  dTJdz  > dTJdr,  one  can  draw  some  interesting  conclusions.  The  regions  above 
and  below  the  center  of  the  laser  spot  have  EiBik.  leading  to  maximum  absorption  there,  if 
the  light  has  s polarization.  (For  p polarization,  maximum  absorption  is  to  the  right  and  left  of 
the  laser  spot.) 


Now  consider  what  happens  if  the  target  is  tilted.  The  geometry  no  longer  has  cylindrical 
symmetry  and  the  magnetic  field  lines  will  be  distorted  into  ellipses,  as  shown  in  Fig.  2b.  No- 
tice now  that  ElBik,  on  the  long  side  of  the  ellipse  for  s polarization,  while  ElBik  on  the 
short  side  of  the  ellipse  for  p polarization.  Thus  absorption  by  self  generated  ion  acoustic  turbu- 
lence favors  s polarization  over  p polarization  if  the  target  is  tilted.  This  may  be  one  possible 
explanation  for  the  fact  that  even  in  experiments  which  show  resonant  absorption,^’  the  ab- 
sorption of  p polarized  light  is  only  about  50%  higher  than  the  corresponding  absorption  for  s 
polarized  light. 

V ABSORPTION  IN  A MAGNETIZED  PLASMA 
QUANTITATIVE  CALCULATIONS 


In  this  section  we  describe  in  more  detail  the  calculations  of  absorption  efficiency  in  a 
magnetized  plasma.  The  steady  state  fiuid  is  described  by  the  four  equations  for  number  densi- 
ty, total  momentum,  electron  energy  and  ion  energy.  They  are 
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Above,  n is  the  electron  number  density,  e}  (EJ  ) is  the  electric  field  squared  of  the  in- 
cident (reflected)  laser  light,  Cfg  is  the  ion  acoustic  wave  induced  energy  exchanged  between 
electrons  and  ions,  v^„  is  the  ion  acoustic  wave  induced  collisional  damping  and  Qg  is  the  elec- 
tron thermal  flux.  These  three  quantities  will  be  specified  shortly. 


The  equations  for  the  spatial  evaluation  of  the  ion  acoustic  waves  are 
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where  d>  is  the  electrostatic  potential  at  wave  number  )c,  V is  the  group  velocity,  which  in  this 
case  is  the  fluid  velocity,  and  0 (*)  is  the  growth  rate^  of  a mode  at  wave  number  k 
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where  5 is  the  magnetic  field.  Equation  (15)  above  is  the  growth  rate  given  in  Ref.  3,  but  with 
ion  Landau  damping  added  on  and  also  with  account  taken  of  the  fact  that  Z ^ 1.  The  mode 
phase  velocity  is  assumed  to  be  in  the  y direction  and  the  gradients  are  assumed  to  be  in  the  x 
direction.  The  configuration  we  adopt  has  the  laser  at  at  •=  — «>  and  the  critical  surface  just  to 
the  right  of  x —0.  Since  energy  flux  is  to  the  right,  the  plasma  temperature  gradient  is  nega- 
tive. Therefore,  as  is  apparent  from  Eq.  (15),  unstable  waves  propagate  in  the  negative  direc- 
tion. 


Coupled  to  the  equations  for  fluid  quantities  and  ion  acoustic  waves  are  the  equations  for 
the  incident  and  reflected  laser  light.  If  9 is  the  target  tilt  angle,  these  equations  are 
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The  remaining  quantities  to  specify  are  Q,  Cy^  and  The  quantity  ^is  given  in  Ref.  3 and 
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where  <u„  is  the  electron  cyclotron  frequency.  The  quantity  Cy^  is  given  in  Ref.  2 for  an  un- 
magnetized plasma.  Using  Ref.  3,  Cy^  can  easily  be  generalized  to  the  case  of  a magnetized 
plasma.  The  result  is 


r r’  1/  3/;)^,mcAt  dT, 

^Te  “ I fl  I £ I 

k ' ' 


2eR|  k\  dx 


ed>  (k) 


(18) 


dn 


Notice  that  the  second  term  in  the  parentheses  above  cools  electrons  since  k,——  < 0.  The 

dx 

explanation  for  this  effect  is  given  in  Ref.  2.  The  one  remaining  thing  to  specify  is  Vg„.  This 
is  calculated  as  in  Refs.  1 and  2,  using  the  theory  of  Dawson  and  Oberman.'*^  However  now 
the  absorption  is  enhanced  because  in  the  configuration  we  examine,  k | E.  The  result  is 
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where  is  the  electron  Debye  length.  The  fluid  system  is  now  completely  specified. 

Equations  (13)-(19)  are  integrated  numerically  back  towards  x — — <»,  subject  to  the 


boundary  conditions  described  in  the  previous  section.  Also 
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“10  ^ at  X -0. 


Typical  solutions  are  shown  in  Figs.  3 a and  b for  a Nd  laser  produced  plasma.  In  both  cases 
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- 20.  « - 10‘  Gauss,  «(x -0)  - 1/2  n„  tf -0.  and  ihe  laser  irradiance  is 
2 X 10*^  W/cm^.  In  Fig.  3a  no  resonant  absorption  is  assumed  to  the  right  of  x -0  and  the  total 
fractional  absorption  is  16%.  In  Fig.  3b,  20%  resonant  absorption  is  assumed,  and  the  fractional 
absorption  is  33%.  In  each  case.  10  values  of  k were  used.  Notice  that  in  each  case  the  ab- 
sorption takes  place  in  a length  of  less  than  10  microns.  Also  the  turbulent  field  strength, 

. 2 

-SL  < ^0  is  less  than  that  calculated  in  Ref.  1. 

Te 

In  Fig.  4 is  shown  a plot  of  flux  limit  / - Q(x  -0)/n»jF^-  (x  -K))  and  fractional  ab- 
sorption as  a function  of  laser  irradiance  between  about  4 x 10*^  W/cm^  and 
1.5  X 10**  W/cm^.  Again  nfx— 0)  — 1/2  , Tj,//",  = 20  , # * 0 , 5 = ( 1 MG  ) 

X (r,  (x -0)/12keV)*/^.  This  functional  form  for  B specifies  a constant  plasma  B as  the 
temperature  is  increased.  Twenty  percent  resonant  absorption  at  the  critical  surface  is  assumed 
here.  Actually  the  flux  limit /should  be  somewhat  higher  than  shown  here,  because  the  elec- 
tron thermal  flux  arising  from  the  20%  assumed  resonant  absorption  was  not  included  in 
Qix  -0).  Figure  4 indicates  the  same  general  result  found  in  Ref.  1,  the  fractional  absorption 
and  flux  limits  are  weakly  varying  functions  of  laser  irradiance,  while  the  electron  temperature 
at  X -0  is  a strongly  varying  function  of  it. 

In  another  series  of  runs,  we  tested  the  dependence  of  fractional  absorption  on  magnetic 
field  strength.  For  /»(x  =0)  - 1/2  0 ‘O.TJT,  - 20,  and  assumed  resonant  absorption  of 

20%  and  a laser  irradiance  of  2 x lo'^  Wlcni^,  the  magnetic  field  was  varied  from  700  KG  to 
2.2  MG.  As  is  shown  in  Fig.  5,  the  relative  absorption  is  about  40%  nearly  independent  of  field 
strength.  However  the  required  length  for  absorption  increases  as  the  field  strength  decreases. 
Thus  if  length  scale  is  an  important  feature  in  limiting  the  absorption,  higher  magnetic  fields 
give  rise  to  higher  absorption. 

In  another  series  of  calculations,  the  density  at  x -0  was  varied  assuming,  «=0, 
/ - 2 X lO'*  Wtcm^,  TJTi  -20,  B - 1 MG,  and  20%  resonant  absorption  at  x -0.  The 
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result  is  shown  in  Fig.  6.  The  results  are  quite  similar  to  earlier  calculations  in  a field  free  plas- 
ma. At  n(x  -H»-0.25  there  is  virtually  no  absorption  (recall  20%  resonant  absorption  was 
assumed).  The  fractional  absorption  however  is  a steeply  rising  function  of  density  at  x -0. 

Notice  that  the  fractional  absorption  rises  to  almost  80%  for  «(x  =0)  “0.7n„. 

In  a final  series  of  calculations  we  examined  the  functional  dependence  of  fractional  ab- 
sorption on  temperature  ratio  for  n(x  M))  =0.5n^,,  I = 2.10*^  Wlcnt^,  B “ IMG  and  also  as-  ] 

\ 

suming  20%  resonant  absorption.  The  results  are  shown  in  Fig.  7.  As  is  apparent,  one  finds  S 

absorption  at  temperature  ratios  as  low  as  5.  One  very  interesting  and  surprising  result  is  the 
maximum  for  low  temperature  ratios.  We  find  a qualitatively  different  behavior  at  low  tem- 
perature ratio  from  that  at  high  temperature  ratio.  At  high  temperature  ratio,  the  turbulent 
spectrum  is  dominated  by  a single  k at  all  x.  For  low  temperature  ratios  however  different  Ac’s 
dominate  the  spectrum  at  different  positions.  This  gives  rise  to  a choppy  behavior  of  n(x)  and 

-— (x)  as  shown  in  Fig.  8. 

‘e 

VI  SCALING  WITH  LASER  WAVELENGTH 

All  of  the  calculations  in  the  previous  section  were  done  assuming  the  laser  wavelength 
was  1.06m.  However  it  is  a simple  matter  to  show  that  these  equations  have  a simple  scaling  so 
that  a solution  for  one  laser  wavelength  can  be  scaled  into  a solution  for  any  laser  wavelength. 

Let  us  assume  that  the  laser  wavelength  increases  by  a factor  of  o,  then  the  equations  in  the 
previous  section  are  all  invariant  under  the  following  scale  transformation 

X — ax 

n — n/a^ 

r—  T 

K—  y 

20 


\ 


(20) 

n — n/a 

B—  B/a 

e<i>  e<i> 

T T ' 

I g t g 

Some  of  the  results  of  this  scaling  transformation  are 

Q-  Qla^ 

^an  *'an/<=‘ 

Cfg  -*  Cfg/a  ^ 

(21) 

y — -y/a 

By  examining  the  equations  for  the  unmagnetized  plasma  in  Ref.  1,  it  is  not  difficult  to 
see  that  these  equations  are  invariant  under  the  same  scale  transformation.  Thus  Figs.  3-8  can 
easily  be  interpreted  for  a CO2  laser  produced  plasma  as  well  as  for  a Nd  laser  produced  plas- 
ma. 


VII  THE  DYNAMICS  OF  THE  INTERACTION 

The  calculations  we  have  performed  in  this  paper  and  in  our  previous  work ' are  all 
steady  state,  that  is  all  d/dt's  are  set  equal  to  zero  in  all  fluid  equations.  As  we  have  discussed 
previously,  one  drawback  with  this  approach  is  that  several  quantities  have  to  be  specified  arbi- 
trarily. One  possible  way  to  specify  these  quantities  in  a physical  way  is  to  solve  for  the 
dynamics  of  the  interaction  rather  than  the  steady  state.  The  way  to  include  anomalous  tran- 
sport and  absorption  into  a dynamic  one  or  two  dimensional  fluid  simulation  is  given  in  Ref.  2 
for  an  unmagnetized  plasma,  and  in  Ref.  3 for  a magnetized  plasma.  By  solving  the  dynamic 
problem,  one  could  then  solve  for  such  quantities  as  0{x  — — «»  ),  n(x  -0 )/«(.,,  etc.  in  terms  of 
boundary  conditions  imposed  far  from  the  critical  surface,  that  is  in  the  solid  target  and  in  the 
vacuum. 
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Fig.  6 — Fractional  abaorption  at  a function  of 
the  lower  shelf  density  n/n 


Fig.  7 — Fractional  absorption  as  a function  of 
dectron  to  ion  temperature  ratio  for  the  case  of 


